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We study the dynamics and thermodynamics of one-dimensional spin-orbital models relevant for 
transition metal oxides. We show that collective spin, orbital, and combined spin-orbital excitations 
with infinite lifetime can exist, if the ground state of both sectors is ferromagnetic. Our main focus 
is the case of effectively ferromagnetic (antiferromagnetic) exchange for the spin (orbital) sector, 
respectively, and we investigate the renormalization of spin excitations via spin-orbital fluctuations 
using a boson-fermion representation. We contrast a mean-field decoupling approach with results 
obtained by treating the spin-orbital coupling perturbatively. Within the latter self-consistent ap- 
proach we find a significant increase of the linewidth and additional structures in the dynamical 
spin structure factor as well as Kohn anomalies in the spin-wave dispersion caused by the scattering 
of spin excitations from orbital fluctuations. Finally, we analyze the specific heat c(T) by compar- 
ing a numerical solution of the model obtained by the density-matrix renormalization group with 
perturbative results. At low temperatures T we find numerically c(T) ~ T pointing to a low-energy 
effective theory with dynamical critical exponent z = 1. 

PACS numbers: 75.10.Pq, 75.30.Et, 05.10.Cc, 05.70.Fh 



I. INTRODUCTION 

In condensed matter systems the coupling between 
different degrees of freedom often plays an important 
role. The electron-phonon coupling, for example, can 
lead to the formation of renormalized quasiparticles, so- 
called polarons^ as well as to phase transitions like the 
Peierls instability^ In recent years, the coupling between 
fcrmionic and bosonic degrees of freedom has also been 
intensely studied in Bose-Fermi (BF) mixtures of ultra- 
cold quantum gases 

Coupled degrees of freedom also seem to be important 
in certain transition metal oxides where the low-lying 
electronic states (termed "orbitals") are not completely 
quenched so that temperature or doping can lead to a sig- 
nificant redistribution of the valence electron density. In 
insulating materials with partly filled degenerate orbitals 
the superexchange between the magnetic degrees of free- 
dom then becomes a function of the orbital occupation. 
This leads to models of coupled spin and orbital degrees 
of freedom as, for example, the Kugcl-Khomskii model 
where the orbitals are represented by a pscudospinr 9 . 
LaMn0 3 ^fl-ia LaTi0 3 ^ii and LaV0 3 or YV0 3f ±&^ 
are well-known examples for compounds believed to be 
described by effective spin-orbital models. They exhibit 
a wide range of fascinating effects ranging from colossal 
magnetoresistancei^ to temperature-induced magnetiza- 
tion reversals ! 18 ' 19 



Common to all these transistion metal oxides is a lift- 
ing of the fivefold degeneracy of the d orbitals into two e g 
orbitals {x 2 — y 2 and 3z 2 — r 2 ) and three t2 g orbitals (xy, 
yz, and xz). This splitting is due to the perovskite struc- 
ture where oxygen ions, O 2- , form octahedra around the 
transition metal ions which are therefore exposed to an 
approximately cubic crystal field. As a consequence, the 
orbitals pointing towards the oxygen ions are energeti- 
cally unfavorable. 

In YVO3 the t2g orbitals are occupied by two electrons 
forming an effective spin 5 = 1 due to large Hund's rule 
coupling. The material is an insulator with an interest- 
ing phase diagram^ -20 ' 22 At temperatures below 77 K 
the system is in a G-type antiferromagnetic (AF) phase, 
i.e., AF in all three directions. In a range of higher tem- 
peratures, 77 K < T < 116 K, the magnetic structure is 
C-type with spins ordering antiferromagnetically in the 
(a, b) plane and ferromagnetically along the c axis. The 
surprising fact that the ferromagnetic (FM) exchange in- 
tegral in this phase is much larger than the AF exchange 
interactions in the (a, b) planed was explained by strong 
orbital fluctuations along the c-axis chains that trigger 
ferromagnetismi^i In the C-type phase a neutron scat- 
tering study revealed that the magnon dispersion along 
the FM c-axis chains consists of two branches. This split- 
ting has been interpreted as due to a periodic modula- 
tion of the FM exchange along these chains caused by 
an entropy gain of fluctuating orbital occupations! 22 ' 23 
Support for an orbital Peierls effect in this material was 
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equation of motion method leading to 
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FIG. 1: (a) Ground state of the FM spin-orbital model, 
Eq. (|l.ip . (b) a spin excitation, and (c) a coupled spin-orbital 
excitation. The two orbitals per site are assumed to be de- 
generate (the splitting is only for clarity of presentation). 



given by numerical investigation o 23 i 24 and a mean-field 
(MF) decoupling approach^ 

However, the dynamics in such systems cannot easily 
be studied numerically and a MF decoupling is unable 
to explain important features of coupled spin-orbital de- 
grees of freedom 2 ^ as can be seen in the following ex- 
ample. Consider the one-dimensional (ID) spin-orbital 
Hamiltonian^— 



H 



Sj+i + x) (tj ■ r j+1 + y) 



(1.1) 



with ferromagnetic superexchange interaction J < 0, 
where Sj and Tj are spin S and pseudospin r operators 
at site j, respectively, and x and y are constants. For 
general x, y the model has an SU(2) <E) SU (2) symmetry 
and exhibits an additional Z2 symmetry, interchanging 
spin and orbital sectors, if x = y. For S = t = 1/2 and 
x = y = 1/4 the symmetry is enlarged to SU(4) — 

In the following we discuss the case S = r = 1/2 
and we choose x and y such that the ground state 
\Fs,Ft) is given by fully polarized spin and orbital sec- 
tors as illustrated in Fig. HJa). Using the equation of 
motion method we find that the state Sj \F S , F T } shown 
in Fig. [IJb) is always an elementary excitation with dis- 
persion cos(q) = |J|(1 +4y)(l — cosg)/4. Analogously, 
the orbital flip is also an elementary excitation with 
u) T (q) = I J|(l+4x)(l — cosg)/4. However, these collective 
excitations are not the only undamped elementary exci- 
tations of the Hamiltonian, Eq. (jl.lj) . In addition, a cou- 
pled spin-orbital excitation S~t~ , as shown in Fig. [TJc) 
may exist. To investigate this issue we again apply the 



Dj{x,y) 



x 4 ) \ S j t j-i + S J T j+D 



(1.3) 



contains terms which lead to a spatial decoherence of 
the excitation. Hence in order to have a fully confined 
spin-orbital excitation, Dj(x, y) has to vanish, which ob- 
viously is the case ifx = y = 1/4. From this we draw 
the conclusion that confined spin-orbital excitations are 
rather the exception than the rule, relying on the par- 
ticular value of the constants. If we introduce the Bloch 
states 



S 3 T J 



\F S ,F T ) 



(1.4) 



the dispersion of the coupled spin-orbital excitation for 
x — y — 1/4 is given by 



^Sr(g) = |J|(1 - COSQ)/2. 



(1.5) 



Thus for x = y = 1/4 we have u>s{q) = oj T {q) = u/s T (q), 
i.e., the dispersions of all three elementary excitations 
are degenerated Interestingly, they all lie within the 
continuum of spin-orbital excitations given by 7(5, p) = 
ws(g/2 + p) + u) T (q/2 — p). The Hamiltonian, however, 
does not allow for a decay of these three elementary ex- 
citations in the ferromagnetic case. For the case of the 
coupled spin-orbital excitation we see from Eq. (jl.3p that 
such a decay becomes possible once we move away from 
the special point x = y = 1/4. Our conclusions partly 
differ from the ones presented in Ref. H3, where the cou- 
pled spin-orbital excitation is considered as a bound state 
below the spin-orbital continuum^ 

There are several ways to generalize the S = r = 1/2 
case to arbitrary spin- and pseudospin quantum numbers. 
If we start again from fully polarized spin and orbital sec- 
tors and only demand that S~t~ \ F$, F t } stays confined, 
we find the condition x = 5(1 — S) and y = r(l — r). 
Another way of generalizing the S = t = 1/2 case to ar- 
bitrary S and r relies on the fact that the Hamiltonian, 
Eq. (|1.1|) . with x = y — 1/4 is equivalent to 



H 



J 



j,3 + V 



(1.6) 



. { with a £ {S, t} is Dirac's exchange operator 



where T> 

for a = 1/2 j2£ A generalization of this exchange operator 
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to arbitrary spin has been discussed by Schrodinger^ 5 - 
For instance for a = 1 the spin (pseudospin) exchange 
operator is given by T>°j 1 = (cTj ■ <x ; ) 2 + a-j ■ eri — 
The Hamiltonian in Eq. (|1.6p with arbitrary spin (pseu- 
dospin) quantum number does not only keep the sin- 
gle spin-orbital flip confined as in the generalization dis- 
cussed above but rather all spin-orbital excitations of the 
type {S~) ms (T~) mT \Fs, F t ) where m s and m T are the 
multiplicities for spin and pseudospin, respectively. Since 
a MF decoupling solution treats the spin-orbital chain 
as two separate chains with effective exchange parame- 
ters determined self-consistently, the physics of coupled 
spin-orbital excitations cannot be captured within this 
approach. 

The purpose of this paper is to study the importance 
of coupled spin-orbital excitations in a spin-orbital model 
with antiferromagnetic superexchange and anisotropic 
orbital exchange. This case is intriguing as spin and or- 
bital degrees of freedom may be expected to be strongly 
entangled^ In fact, it has been shown that compos- 
ite spin-orbital excitations have to be analyzed together 
with spin waves in systems with active e g orbitals, such 
as for instance KCuF3i 37 i 38 This follows from the non- 
conservation of the orbital flavor in hopping processes 
which implies that spin excitations are not independent 
and may occur in general together with an orbital flip. 
Here we will consider an anisotropic generalization of the 
spin-orbital model with parameters x, y such that 
the spins still order ferromagnetically in the ground state. 
The orbital sector, however, will no longer be in a fully 
polarized state due to the AF superexchange which fa- 
vors orbital alternation. Independent spin, orbital, and 
coupled spin-orbital excitations of collective type, as dis- 
cussed above, therefore can no longer exist. We will fo- 
cus, in particular, on the question how spin excitations 
are modified by the presence of orbitals in this case. 

The paper is organized as follows: In Sec.|H]we present 
a generalization of the spin-orbital model, Eq. (jl.ip . to a 
model with anisotropic orbital exchange. For the extreme 
quantum limit of the orbital sector interacting via an XY- 
type coupling we then derive an effective BF model which 
resembles models considered in the context of ultracold 
BF gases. By using a density matrix renormalization 
group algorithm applied to transfer matrices (TMRG) we 
cxemplarily investigate numerically the crossover from 
AF to FM correlations. In Sec. [HI] we discuss the MF 
decoupling approach. We allow for a dimerization in 
both sectors and discuss the obtained MF phase diagram. 
In Sec. IIVI we summarize the results for the dynamical 
spin structure factor S(q, ui) obtained within the modi- 
fied spin-wave theory (MSWT) for the uniform FM spin 
chain] 39 i 40 In Sec. |V] we formulate an approach where 
the coupling between spins and orbitals is treated per- 
turbatively. The approach is based on representing the 
spins by bosons using the MSW T 41 ' 42 and the orbitals by 
Jordan- Wigner fermions. Finally, in Sec. lVIl we consider 
the effects of coupled spin-orbital degrees of freedom on 
the thermodynamics of the system. We focus, in partic- 



ular, on the specific heat as a function of temperature 
and compare perturbative results with numerical data 
obtained by TMRG. In Sec. IVIII we summarize and dis- 
cuss our results. The Appendix provides details of the 
perturbative approach. 



II. SPIN-ORBITAL MODEL AND MAPPING 
ONTO A BOSON-FERMION MODEL 

A. One-dimensional spin-orbital model 

We will focus here on the physical situation realized 
in the vanadium perovskites, such as YVO3, where the 
superexchange interactions are antiferromagnetic. In 
YVO3 the two d electrons occupy the lower lying ti g or- 
bitals while the e g orbitals are empty. From electronic 
structure calculations^ 3 - - — it is concluded that the t-2 g 
orbitals are split into a lower-lying xy orbital level and 
a higher-lying doublet of xz and yz orbitals. Therefore 
the xy orbital will always be occupied by one electron, 
controlling the AF correlations in the (a, b) planes. The 
large Hund's coupling Jh (normalized to the interatomic 
Coloumb interaction U) present in YVO3 will support 
parallel alignment of electronic spins at V 3+ ions in d 
configurations^ leading to S = 1 spins. Therefore, the 
remaining electron will be placed in one of the two other 
orbitals {xz,yz} which constitutes the t — 1/2 orbital 
degree of freedom. On this basis, an effective spin-orbital 
superexchange model for YVO3 with spins S = 1 has 
been derived^ Here we shall study the ID spin-orbital 
model extracted from it for the c-axis.— 

The simplest Hamiltonian for the c-axis FM chains in 
YVO3 taking Jh into account is given by Eq. with 
J > 0, x = 1, and y = 1/4 — 711. Here 7h is proportional 
to Hund's coupling Jh, supporting FM correlations in 
the spin sector. For 5 = 1 and realistic values of Hund's 
coupling for vanadates, 7h ~ 0.1, numerical investiga- 
tions of this model showed strong but short-ranged dimcr 
correlations in a certain finite temperature range caused 
by the related entropy gain although the ground state is 
uniformly FM<2i The same Hamiltonian was also studied 
using a MF decoupling schemed Within this approach a 
finite temperature phase with dimer order in both sectors 
was found. However, as discussed in the introduction, the 
MF decoupling approach has severe limitations as it does 
not take the coupled spin-orbital dynamics into account. 

We start from a generalization of Eq. (II. ip which reads 

H S T(r) = jJ^(S j -S j+1 +x)([T j -T j+1 ] T +y) , (2.1) 

3 

with 

[Tj ■ r j+ i] r = Tj ■ r j+1 - Tt*t* +1 . (2.2) 

The calculations presented here are valid for general S 
and x, but we will, unless stated otherwise, only address 
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FIG. 2: (Color online) Nearest-neighbor spin and orbital cor- 
relation functions for the spin-orbital model (|2.1|) with r = 1 
as a function of temperature T in units of J (we set &b = 1). 
In both panels y = 0.4, 0.3, 0.2, 0.15, 0.1, 0.05, 0.0, -0.1, 
—0.2, —0.3 in arrow direction. The spin correlations switch 
from AF to FM at y = 0.1 (dashed lines). The dotted lines 
in the upper (lower) panel correspond to the limiting values 
1 and —1.4015 ( — and l/7r), respectively. 



the case S = 1, x = 1 relevant for YVO3 in the following. 
For r = 1 the pseudospin sector reduces to an XY model. 

In Fig. [3 the nearest-neigbor spin and orbital corre- 
lation functions for T = 1 as a function of temperature 
for various parameters y obtained by TMRG are shown. 
This method allows us to obtain thermodynamic quan- 
tities for ID quantum systems directly in the thermody- 
namic limit For y < 0.1 the ground state has ferro- 
magnctically aligned spins. The phase transition between 
the fully polarized FM state and a state with AF spin cor- 
relations at y w 0.1 is first order. In the limit y 3> 1 the 
value (SjSj+i) ~ —1.4015 for a Haldane 5 = 1 Heisen- 
berg chain is reached^ while the orbital correlations ap- 
proach (TjTj +1 + TjTj +1 ) —> 1/tt. We note that the FM 
ground state is lost at y ~ 0.1 both in the model with 
r = 1 investigated here as well as in the model with an 
isotropic pseudospin sector (r = 0) i 24 i 51 For the model 
with r = 1, however, we have a direct phase transition 
from the FM to the Haldane phase while for the isotropic 
model an orbital valence bond phase is intervening be- 
tween these two phases^ 

The isotropic model (|2.ip . T = 0, has also been in- 
tensely studied for S = r = 1/2. Here the phase diagram 
is more complex than in the S — 1, r = 1/2 casei 52 ' 53 



For x = 1 the FM spin state is again found to be stable 
for y < 0.1. However, now the transition at y ~ 0.1 is to 
a gapless "renormalizcd SU(4)" phase followed by a fur- 
ther phase transition at larger y into a dimer phase. The 
phase with ferromagnetically polarized orbitals is absent 
because (Sj ■ Sj+i +x) > for x = 1 but is again present 
for large y if x < In 2 — 1/4. 



B. Boson- fermion model 

For the ID spin-orbital model, Eq. (|2.ip . at the point 
r = 1, we will now derive an effective BF model which 
will be used as a starting point for the perturbative 
approach. First, applying the Jordan- Wigner transfor- 
mation the orbital part is mapped onto a free fermion 
model (for r ^ 1 the pseudospins map onto interacting 
fcrmions). The spin part of the spin-orbital Hamiltonian 
(|2.1|) will be represented by bosons. Concentrating on 
the case where the spin part is ferromagnetically polar- 
ized in the ground state, we can treat the spin sector by 
the MSWTi 41 i 42 To this end, we introduce bosonic op- 
erators by a Dyson-Maleev transformation. If we retain 
bosonic operators only up to quadratic order we end up 
with 

H = 'Hsr(l)-JN{S 2 + x)y^H + H 1 . (2.3) 
Here Hq is already diagonal 

H =Y / "B(k)blb k + X>(9)4/ 9 , (2.4) 



with /j and f q {b\ and bk) being the fermionic (bosonic) 
creation and annihilation operators, respectively. The 
magnon dispersion is given by 



wbO) =2J5|y|(l-cosfc), 
and the fermion dispersion reads 

w f(?) = J(S 2 + x) cos q . 



(2.5) 



(2.6) 



The spinless fcrmions fill up the Fermi sea between the 
Fermi points at fcj? = ±7r/2. 

For FM spin chains usual spin-wave theory has to be 
modified by a Lagrange multiplier \i acting as a chemical 
potential which enforces the Mermin- Wagner theorem of 
vanishing magnetization at finite temperature 4 ^ 



,5' 



N 4-^ 



(2.7) 



Thermodynamic quantities calculated with this method 
are in excellent agreement with the exact Bethe ansatz 
solution for the uniform chain as well as with numerical 
TMRG data for the dimcrized FM chain for temperatures 
up to T ~ I J c ff IS* 2 , with J c ff being the effective exchange 
constant of the model under consideration ! 26 ' 41 ! 42,54 
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The interacting part couples bosons and fermions and 
reads 

Hl = jf E "Mki,k2,q)blb k j}f kl _ k2+q , (2.8) 

ki ,k 2 ,q 

with the vertex 

WBF(fci, k 2 , q) = JS [cos(/c 2 - q) + cos(fci + g) 



cos(fci — /c2 + q) — cos g] 



(2.9) 



The Hamiltonian H = Hq + Hi with Ho and Hi given 
by Eqs. (|2.4p and (|2.8|) . supplemented by the constraint 
(|2.7p , is an effective BF representation valid at low tem- 
peratures. We will investigate this model in Sec. [V] treat- 
ing the BF coupling perturbatively. 



III. MEAN-FIELD DECOUPLING 

The spin-orbital model (|2.1j) contains rich and inter- 
esting physics. A first attempt to understand the prop- 
erties of the model is to apply a MF decoupling which 
neglects the coupled spin-orbital degrees of freedom and 
treats the spin-orbital chain as two separate chains with 
effective coupling constants which have to be determined 
self-consistently. Note, however, that this treatment does 
not involve site variables as in the classical Wciss-MF 
theory but takes the correlations on a bond as relevant 
variables. Interestingly, these expectation values never 
vanish, which makes them useful particularly in cases 
without long-range order. 



A. Decoupling into spin and orbital chain 

Applying a MF decoupling and allowing for a dimcr- 
ization in both sector o 26 i 27 we obtain from Eq. (|2.1[) 

n S r(r)~Hf F + H™ F (r), (3.1) 

with the spin and orbital Hamiltonians 

N 



Hf F = J 3 J2{ l + (- l Y 6 s}S r S J+ i, 



N 



(3.2) 



i? r MF (r) = J T J2{1 + (-l) j Sr} [Tj ■ T i+1 ] I 

3=1 



Within this approximation the effective superexchange 
constants and dimerization parameters are given by 



Jr = J ■ 



St = 



A 



ss 



A+ + 2y A 
Js = J 7, ' d s 



A ^ + 2X (3.3) 



where we have defined 
- IS** ■ S 



»2j 

'V 



2j+l) ± (S23 ■ Sij-i) , 
T 2] + l] r ) ±([T 2 j -T 2j-l] T ) ■ 



(3.4) 



Here A~ a with a = S (a = t) is an order parameter for 
the spin (orbital) dimerization, respectively. Thus, the 
exchange constants and dimerization parameters for each 
sector are determined by the nearest-neighbor correlation 
functions in the other sector, making a self-consistent 
calculation necessary. In the following we want to solve 
Eqs. (|3.2p - (|3.4p for the spin exchange being effectively 
FM, i.e. J s < 0. 



B. Dimerized orbital correlations 

Numerical investigations of the model with isotropic 
orbital exchange, Hst(0), have shown orbital-singlet for- 
mation in the ground stated for y > 0.1. Moreover, 
although the ground state consists of a fully spin polar- 
ized FM state for y < 0.1 with AF orbital correlations, 
it has been shown that a tendency towards orbital sin- 
glet formation is still present but has to be activated by 
thermal fluctuations^ In Ref. [26] the model (|3.ip was 
studied in the FM regime with x = 1, y = \ — 7h, 
T = and 7h = 0.1 in order to address the ques- 
tion whether this orbital-Peierls effect can be captured 
within a MF decoupling approach. A dimerized phase 
for 0.10 < T/J < 0.49 (we set k B = h = 1) was found 
with the dimerization amplitude in the spin sector being 
much larger than in the orbital sector. 

We now want to compare this result with the case 
where we set T = 1 in Eq. (|3.1[) so that the self- 
consistent Eqs. (|3.3p can be solved analytically by ap- 
plying a Jordan- Wigncr transformation and MSWT. In- 
troducing fermionic operators /■ ^ if the index j is even 

and J~j o ^ 3 1S °°-d for the pseudospins, we rewrite 
^jMF = H^ F (T = 1) in Fourier representation. Finally 
introducing new fermionic operators 4>^ and (p^ which 



diagonalize the Hamiltonian H^r , we find 

i? r MF = E-F F (9,^)(^ + ^ 9 )' 



with the fermionic dispersior 
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(3.5) 



^ (9, S T ) = J T J cos 2 q + SI sin 2 q . (3.6) 



We can now calculate A^ r , as given in Eq. 
straightforwardly and obtain 

A _ _ 2Sr^{2n F [ojf F (q,S T )}-l}sm 2 q 



N 



cos 2 q + S 2 . sin q 



A+ T + 2y 



2 v {2n F [< F ( g ,^)]-l}cos 2 g 

I 2 , „ • 2 

cos q + o£ sin q 



(3.7) 



where n-p(x) = {cxp(/3x) + 1} 1 is the Fermi function 
and = 1/T. 



C. Dimerized spin correlations 

Next we turn to the spin part of Eq. (|3.ip to which 
we apply the MSWT i 41 ' 42 We introduce two bosonic op- 
erators tfjl [bj'l] for j even [odd] by means of a Dyson- 
Maleev transformation. Retaining only terms bilinear in 
the bosonic operators we can diagonalize the resulting 
Hamiltonian by a Bogoliubov transformation leading to 

Hf F = £ {u%?_ (k, 6 s )a{a k + u%? + (k, S s )^ k } 

k 

+ JsNS 2 , 

(3.8) 

with the two magnon branches 



u%X(k, 5 S ) =2\J S \S ( 1 ± ^cos 2 fc + 5 2 sin 2 k) . (3.9) 



The constraint of vanishing magnetization at finite tem- 
perature (|2.7|1 now reads 



S = ^£MC(Ms)]+n B [C + (Ms)]}, (3.10) 



where ns(x) = {cxp{/3x'} — 1} 1 is the Bose function and 



C ± (k,S s ) = ^l(k,6s)- l x(S s ). 



To calculate the nearest-neighbor correlation functions 
B± = (Sj ■ Sj±i\ it is necessary to go beyond linear spin- 
wave theory. Taking terms of quartic order into account 
and using Eq. (|3.10|) we obtain 2 ^ 

£±=(^£/±(Ms) £ an B [C B (k,S s )}\ 



Here we have defined 



/±(Ms 



cos 2 k ± 5s sin k 



cos 2 k + 5 S sin 2 k 



(3.11) 
(3.12) 



From these expressions we can obtain A ss which, com- 
bined with Eq. ([3TT]) . allows us to solve Eqs. (pTTj) - (pHl) 
self-consistently. 



D. Mean-field phase diagram 

We first discuss the ground state phase diagram of the 
Hamiltonian (|3.1[) for V = 1. Depending on the sign of 
the effective coupling constant Js we find (SjSj+i) = 
1, —1.4015 with the latter value being the approximate 
result for the 5 = 1 AF Haldane chain. In the following, 
we restrict our discussion to — 1 < x < 1.4015 so that 



0.4 - (a) T i^ 



- tricri tical 
point 



uniform 




0.3 

T/J T/J 



T 2 /J 



FIG. 3: (a) Phase diagram of the Hamiltonian H3.1[) with 
r = 1 and x = 1 in mean-field decoupling. The shaded area 
represents the dimerized phase. The phase transition at T2 
is first order whereas the transition at T\ is of second order. 
The two transition lines merge at the tricritical point y tp . (b) 
Dimerization parameters 5s and S T for x = 1 and y — 0.14. 
The lines are guides to the eye. The shaded area marks the 
temperature range where the dimerization is nonzero. 



(SjSj+i) and J T = J((SjSj+i) + x) always have the 
same sign. For the orbital sector we obtain, on the other 
hand, (tjTj +1 +t^t^ +1 ) = ±l/ir. y > 1 /it implies J s > 
and the ground state is therefore certainly AF (Haldane 
phase) whereas Js < for y < —1/tt leading to a FM 
state. In the regime —l/ir < y < 1/tt the self-consistent 
equations have two solutions with energies Eq F w (1/tt + 
y)(-lA015 + x) and^ M = (-l/n + y)(l+x) and a first 
order phase transition between the FM and AF states 
occurs where the energies cross. For the case x = 1 we 
are focussing on here, this happens at y c w 0.212 and the 
FM state is stable for y < y c . Compared to the numerical 
solution where y c rj 0.1 (see Fig. [5]) the range of stability 
of the FM state is therefore increased in the MF solution. 

Next, we investigate the possibility of a finite tem- 
perature dimerization for x = 1 in that part of the 
phase diagram where the ground state is FM. As shown 
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in Fig. (3Ja) we find that a dimerized phase at finite 
temperatures does indeed exist in MF decoupling for 
y tp ~ 0.128 < y < y c « 0.212 where y tp denotes the 
tricrictal point. As in the model with an isotropic pseu- 
dospin sector^ the temperature range where the dimer- 
ized phase is stable depends on y. At the onset temper- 
ature T\ the phase transition is of second order whereas 
at the reentrance temperature T2 it is of first order, see 
Fig. EIb). As in the case T — 0, the dimerization in 
the spin sector is always much larger than in the orbital 
sector. 

As pointed out before, the MF decoupling suffers from 
severe limitations and it is expected to be an even worse 
approximation in the extreme quantum case T = 1 than 
in the case T = studied previously^ In particular the 
coupling between spin and orbital degrees of freedom is 
completely lost within this approach. In the following 
sections we will therefore develop an alternative pcrtur- 
bative treatment of the spin-orbital coupling. 



IV. DYNAMICAL SPIN STRUCTURE FACTOR 
FOR THE UNIFORM FERROMAGNETIC CHAIN 

In order to investigate coupled spin-orbital degrees of 
freedom and, in particular, their implications on the spin 
dynamics of the spin-orbital chain, a detailed under- 
standing of the spin dynamics of a FM chain is useful. 
We shall avoid the complications of the dimerized chain 
and focus our study on the uniform ID ferromagneti^ 
In doing so we neglect the coupling between spin and 
pseudospin operators for a moment and consider 
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with Js < 0. It is well-known that MSWT does not 
respect the SU(2) symmetry of the FM Heisenberg chain 
Eq. (|4.1[) . We therefore directly calculate the full spin 
correlation functio n 39 ' 40 

G(r,T) = -(T[Sj(P)-S j+ r(r)]) . (4.2) 

In Fourier space we obtain 

1 ^ 1 - e"' 36 "^' 

G(q,u v , B )= -J2(l +riB [{(k)})n B [aq-k)}- ttt 

N ^ - e 9 (fc) 

(4.3) 

where we have used the bosonic Matsubara frequencies 
w^b and e q (k) = £(fc) — £(q — k) with k G [— 7r,7r]. The 
reduced magnon dispersion reads 

C(fc) = 2 J S S(1 - cos k) - (j,. (4.4) 

In Fig. IDJa) the dynamical spin structure factor, 

S(q, u) = 2n B (-w) ImG rGt (<?, w), (4.5) 



FIG. 4: (Color online) (a) Dynamical spin structure factor 
S(q,u) as obtained for T/\J S \ = 0.1 and < q < tt. The 
dashed line indicates the upper boundary of the two magnon 
continuum. The dots are projections of the peak positions 
onto the (q, uu) plane. They are connected by the dotted line 
which is a guide to the eye. (b) Dynamical spin structure 
factor S(q,co) for the same parameters at q = 4ir/5 (solid 
line) and the corresponding density of states (dashed line). 

is shown for the uniform FM chain at uj > 0, where 

ImG ret ( g , W ) =^£( X + n B [C(fc)])n B [C(3 - *0] 

k (4.6) 
x (e-P e «W - l) 5{cj - e q {k)) 

is the imaginary part of the retarded Green's function 
obtained from Eq. (|4.3p by analytical continuation. Up 
to a factor of 27r, as a matter of definition, we obtain 
the result previously given by Takahashi.— The structure 
factor fulfills detailed balance, S(q,u>) = S(q, —uj). 



The symbols in Fig. 0] show the peak positions projected 
onto the (q, ui) plane. They follow the reduced disper- 
sion Eq. (|4.4|) . Also shown in Fig. [4ja) as a dashed 
curve is w™ ax = 4| J$ \S sin | corresponding to the up- 
per boundary of the two magnon continuum e q (k) above 
which S{q,u) is zero in this approximation. 

At the edge of the two magnon continuum S(q, u>) has 
a singularity. In Fig. U£b) the dynamical spin structure 
factor for the same parameters as used in Fig. SJa) is 
shown at q = Att/5 together with the density of states 

which is given by p q (u>) = 1/ ^/(w™ ax ) 2 — uj 2 . Right be- 
low the singularity at w™ ax the density of states to lowest 

order reads p g (w™ ax — Suj) ~ 1/VSlo, i.e., S(q,tu) shows 
a square root divergence at the upper threshold. If the 
edge singularity and the central peak arc well separated 
then the spectral weight of the edge singularity is much 
smaller than the spectral weight of the central peak. If, 
on the other hand, the edge singularity is close to the 
central peak then the shape of the latter is strongly af- 
fected by the occurence of the edge singularity. In this 
case the edge singularity gives a significant contribution. 

It is instructive to analyze S(q, ui) in the limit of small 
q. If the edge singularity and the peak of the struc- 
ture factor are well separated, the lineshape of the peak 
can be obtained approximately. To this end, for small 
q but I Js\S 2 q/T 3> 1 we only retain the leading terms 
of Eq. (|4.3[) . Performing a saddle point approximation 
to lowest order we find S{q,oS) ~ tib(— u))(a(q, u>) — 
a(q, —oj)), with 

\Js\Sq 

a(q, w) « 2S * . (4.7) 

(w - J s Sq 2 ) 2 + (^fa) 

This Lorentzian lineshape is only valid for low temper- 
atures. Here £ « \Js\S 2 /T is the correlation length in 
the low-temperature limitj^r— Finally, we want to stress 
that for T — > the peaks will reduce to 5-functions, i.e., 
only thermal broadening is included in this approxima- 
tion. 



V. PERTURBATION THEORY 
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FIG. 5: (Color online) Magnon dispersion u)b(?) (dashed line) 
for y = 1 and fermionic two-particle continuum (shaded area) . 



presence of the fcrmions for q < 7r/2, since it can not cou- 
ple to these degrees of freedom. However for higher mo- 
menta the spin wave may couple to the fermionic degrees 
of freedom and thus a broadening of S(q, u>) should occur. 
Moreover by choosing different values for y the point at 
which the magnon enters the fermionic two-particle con- 
tinuum is changed. Thus the momentum at which the 
spin wave is affected by the coupling to the fermionic de- 
grees of freedom depends directly on the parameter y. 

These arguments give the qualitatively correct picture, 
i.e., we find indeed that the coupling of the magnon to 
the fermionic degrees of freedom has strong effects on 
the dynamical spin structure factor at intermediate and 
high momenta and that the onset of these effects can be 
well estimated by our simple argument. However, there 
are also certain aspects which can not be captured within 
this picture. For instance, for 2S\y\ > (S 2 + x) it suggests 
that the spin wave may leave the fermionic two-particle 
continuum at a certain momentum qi and thus should be 
unaffected by the BF coupling for q > qi. The detailed 
calculation, however, reveals that this is not true because 
the spin wave decays into a fermionic particle-hole and a 
remaining spin wave as will become clear in the following. 



In this section we intent to go beyond the MF decou- 
pling approach treating the influence of the BF inter- 
action, Eq (|2.8|) . on the spin- wave dispersion pcrturba- 
tively. Naively one would expect that the magnon should 
be able to couple to the fermionic degrees of freedom if 
it lies inside the fermionic two-particle continuum. The 
upper and lower boundary of the latter are given by 

ei? ax (<z) - 2J(S 2 + x)sm(q/2), 

ef n (q) = J(S 2 +x)smq, (5.1) 

respectively. The continuum and the magnon dispersion 
cjb(q) for y = — 1 are shown in Fig. [5] One would there- 
fore expect that in this case wb(<?) is unaffected by the 



A. General formulation 

Here we want to study the Hamiltonian H given by 
Eq. (|2.3[) . with its nonintcracting part Hq and interact- 
ing part Hi defined by Eqs. (|2.4|) and (|2.8[) . by treat- 
ing the BF interaction perturbatively, i.e., in the limit 
|ar|, \y\ S> 1. As explained in the appendix, we start by 
performing a MF decoupling for the interaction H\ . Cor- 
rections to this solution are then taken into account per- 
turbatively. Here we adress the bosonic Green's function 
at zero temperature 

S B (q,t) = -i(T t [b q (t)bl(0)]). (5.2) 
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FIG. 6: Diagrams which contribute to a renormalization of the magnon in a perturbation theory: (a,b) Diagrams with momen- 
tum exchange between the magnon and the fermions, and (c) diagram without momentum exchange. All diagrams are second 
order. Fermionic propagators are shown by solid lines, whereas bosonic propagators are shown as dashed lines. 



ki + k 2 




FIG. 7: Second order diagram for a system of interacting 
fermions with momentum exchange. 



In Fig. [6] all distinct, connected diagrams beyond the MF 
decoupling up to second order are shown. 

We calculate the Green's function from the Dyson 
equation 



with 



C(«), 



(5.3) 



(5.4) 



where £(q) is the reduced magnon dispersion defined 
in Eq. (jt?4|) with J s = J(y - and the self- 

energy S(g, u>) is approximated by the proper self-energy 
£2(9,0;) obtained by summing up the diagrams which 
can be composed of the diagrams shown in Fig. [5] 

The diagrams shown in Figs.[6]Ja) and[6]^b) are of par- 
ticular interest because they are the lowest order dia- 
grams where bosons and fermions exchange momentum. 
They describe the part of spin-orbital dynamics which 
cannot be captured within the MF decoupling approach 
discussed in section Hill The diagram shown in Fig. [SJb) 
has to be thermally activated, i.e., it does not give any 
contribution at zero temperature. The same is true for 
the diagram shown in Fig. [6]^c). Thus at T = the only 
second order diagram which contributes to the self energy 



is the one shown in Fig. [S^a) leading to 

pj2 Yl ^BfO?^!,^) 



s bf(<?^) 



e[uj F {q - h + k 2 )}e[-u} F (k 2 )] 



(5.5) 



w-fi+(fci,fc 2 ) + «0+ 
where we have abbreviated^ 

n±(fci, k 2 ) EE ±C(fci) + cj F (q - fcl + k 2 ) - UJ F (k 2 ) , (5.6) 

with ki,k 2 £ [— 7r, 7r]. 

For systems of interacting fermions we know that per- 
turbation theory in one dimension often leads to infrared 
divergencies 1 58 ' 59 Such divergencies occur, for example, 
for the fermionic analogon of the diagram with momen- 
tum exchange, see Fig. [7] These problems can be over- 
come by the Dzyaloshinski-Larkin solution or bosoniza- 
tion techniques. For the model considered here, however, 
we find no divergencies within the considered diagrams. 
One reason for this behavior is a lack of nesting. While 
for a fermionic interaction as shown in Fig. [JJ all the 
dispersions in the denominator of Eq. (|5.5p arc approxi- 
mately linear at low energies here one of the dispersions 
is approximately quadratic so that nesting only occurs 
for singular points. As a further check, we have evalu- 
ated the integrals in Eq. (|5.5p for a constant vertex at 
small q and uj and did not find any infrared divergencies. 

For finite temperatures the Matsubara formalism can 
be applied straightforwardly The self-energy Eq. (|5.5|) 
now reads 



w| F (g, h,k 2 ) 



1 



N 2 g^iu^B-SltfaM) 



xA/j^ B (fci, k 2 , Uu,B,T)N' FiF (q, ki,k 2 ,u v ,B,T), 

(5.7) 



where we have abbreviated 

A^ B (fci, fc 2 ,uv B ,T) 
Af F , F (q, ki,k 2 ,oj v ,B,T) 



«B[C(fci)] ± n F [uj F (k 2 )}, 
n F [uj F (q -ki + k 2 )] 
-n F [wF(*a) - C(fci)]. (5-8) 



i n 




FIG. 8: (Color online) Perturbative results for the BF model at zero temperature with x = 1. In the left (right) panel y — — 1 
(y = —2), respectively. (a),(c) S(q, lj) with the inset showing the region for which ImEg F (g,u) is nonzero as a shaded area 
and the renormalized spin-wave dispersion as a dotted line. While S(q, ui) is sharply peaked at low momenta, a significant 
broadening occurs at higher q. Moreover we find additional structures which, as explained in the text, are due to coupled 
spin-orbital excitations. (b),(d) — ImSg F (g,tj) as given in Eq. (|5.5p for the corresponding values of q shown in (a) and (c), 
respectively. The coupled spin-orbital excitations show up as peaks and edges in ImEg F (g,o;) (notice the logarithmic scale). 



Note that at finite temperatures both, the reduced spin- 
wave dispersion as well as the fermionic dispersion 
is renormalized due to the MF decoupling applied to 
Eqs. (|2.8[) . The respective expressions are given in the 
appendix, see Eqs. (jA.ip and (| A.2|) . 

At finite temperatures also the diagram shown in 
Fig. EJb) contributes and is given by 



x (1 +Afp B {ki,k2,u;„ i B,T))J\f FtF (q,k 1 ,k2,Uv,B,T). 

(5-9) 



B. Dynamical spin structure factor 

Below we present the results obtained by summing up 
the diagrams shown in Fig. [5] in a Dyson series, but re- 
placing the external legs by the SU(2) symmetric function 
given in Eq. (|4.3|) . While the perturbative results can, 
strictly speaking, only be valid for \x\, \y\ ^> 1 we extend 
the results here to more physical values |a;|, \y\ ~ C(l) 
where we still expect perturbation theory to give at least 
a qualitatively correct picture. Numerical results ob- 
tained for the dynamical spin structure factor within this 
perturbative approach are shown in Fig. [5] for T = and 
y = —1 and y = —2. In both cases S(q,uj) is sharply 
peaked at small momenta whereas a significant broad- 
ening occurs at higher momenta. Note, that within the 
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FIG. 9: (Color online) Renormalized magnon dispersions 0J q 
in the FM chain for x = 1 and selected values of y. Inset: 
The most pronounced Kohn anomalies occur at q near n/2. 
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FIG. 10: (Color online) Magnon linewidth F q (FWHM of 
S(q,u)) at zero temperature (data points), as obtained for 
x = 1 and the representative values of y indicated in the plot. 
The lines are guides to the eye. 



MSWT S(q,cu) is always a <5-function for the pure spin 
model at T — 0, i.e., the broadening here is solely due to 
the coupling to orbital excitations. 

By extracting the central peaks of the dynamical spin 
structure factor at various momenta, we obtain the renor- 
malized spin-wave dispersion uj q within the perturbative 
approach. The result of this is shown in the insets of 
Fig.|S{a,c) and in more detail in Fig.[5J The magnon dis- 
persion is renormalized and small kinks are visible close 
to q = 7r/2, which may be interpreted as Kohn anoma- 
lies (see below). The inset of Fig. [5] shows the Kohn 
anomalies with a higher resolution. For itinerant ferro- 
magncts Kohn anomalies are well-known. Here the inter- 
action between the spins of localized ions is mediated by 
an exchange with the conduction electrons <£2r— These 
Kohn anomalies can thus be used to gain information 
about the Fermi surface of the conduction electrons i^2r— 
However to the best of our knowledge Kohn anomalies in 
the spin- wave dispersion for insulating materials have not 
been adressed so far. As we will show below, the Kohn 
anomalies in our case are caused by coupled spin-orbital 
degrees of freedom. 

Apart from extracting the effective spin-wave disper- 
sion from S(q,oj), we also want to discuss the magnon 
bandwidth (full width at half maximum (FWHM)) T q 
of the central peaks. A broadening of the zero temper- 
ature peaks occurs whenever the imaginary part of the 
self-energy, 

ImS^few) = --^j w bf(<7, k 1 ,k 2 )Q[-u} F (k 2 )] 

x <d[oj F (q-k 1 + k 2 )}5(Lj-n+(k u k 2 )), 

(5.10) 

is non-zero at u> = Q^(ki,k 2 ). The contributions within 
the sums are now determined by the argument of the 5- 
function as well as by the constraints given by the Heav- 
isidc functions. This procedure, for a given set of param- 



eters x, y, and S, effectively yields a region within which 
the spin wave may scatter on fermion pairs. Henceforth 
we call this region the BF continuum. The BF continuum 
is shown in the insets of Figs.jS^a) and[8lc) as shaded ar- 
eas. The upper boundary of the BF continuum (solid 
lines), which is periodic with a period of 2tt, is given by 
-4JS(y - 1/tt) + 2J(S 2 + x) for q = 0, and decreases 
monotonously from this value with increasing \q\. The 
lower boundary (dashed lines) is periodic with a period 

Of IT. 

To obtain the FWHM of the structure factor, the mag- 
nitude of the contributions to the sums given in Eq. 
(|5.10p are essential. Here not only the (hi, ^-region 
which contributes to the summation but also the magni- 
tude of the vertex wbf is of importance. We observe that 
the vertex is small at small momenta but increases at in- 
termediate and high momenta. This leads to a strong 
increase of the magnitude of the imaginary part of the 
self-energy as shown in panels (b) and (d) of Fig. [5J From 
the insets of Fig. [8] it becomes clear that the spin-wave 
dispersion enters the BF continuum depending on y. For 
higher values of \y\ the spin wave enters at lower mo- 
menta. However, since the vertex gives smaller contribu- 
tions at smaller momenta, the broadening of the central 
peaks of the dynamical spin structure factor turns out 
to be smaller the smaller the momenta are at which the 
spin wave enters the BF continuum. This can be seen in 
Fig. [10] where the magnon linewidth T q is shown. The 
onset of a finite T q signals the entrance of the spin-wave 
dispersion into the BF continuum and depending on the 
momentum at which the entrance occurs the increase of 
T q is either smooth (entrance at low momentum) or steep 
(entrance at high momentum). In addition, we observe 
that T q has a maximum at the boundary of the Brillouin 
zone for stronger interactions (y = —0.75 and y = — 1 
in Fig. 1101 respectively), whereas for smaller interactions 
we observe the maximum at smaller momenta followed 
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FIG. 11: (Color online) Effective dispersion of the spin 
wave uj q for y = — 1 (red solid line) together with coupled 
spin-orbital excitations deduced from —ImT,2(q,uj), shown 
by (green) dots within the Bose-Fermi continuum (shadded 
area) . 



by a decrease of the FWHM towards the zone boundary 
(y = -2 in Fig.PTDj). 

Interestingly, the coupling to the orbital degrees of 
freedom does not only give rise to a featureless broaden- 
ing of S(q, uS) but produces additional structures. These 
additional structures are most obvious in Fig.[5Ja). From 
Fig. [Sib) it becomes clear that these structures are dom- 
inated by local extrema as well as edges in the imagi- 
nary part of the self-energy. Eq. (|5.10p shows that such 
extrema can occur if fi+(fci,fc 2 ), Eq. (|5.6p . becomes sta- 
tionary as a function of the momenta k\ , hi as long as the 
Heaviside functions in Eq. (|5.10p for these momenta are 
non-zero. The position of the local maxima in the imag- 
inary part of the self-energy is therefore approximately 
given by the values fl^(k\, fo) at these stationary points. 
These values correspond to the energy of spin-orbital ex- 
citations into which the initial spin wave can decay, see 
Fig. [|2a) and which are stable against small redistribu- 
tions of momenta. 

We conclude that while we do not have completely 
sharp spin-orbital excitations any more as in the Hamil- 
tonian (|1.6p with FM exchange considered in the intro- 
duction, there are still characteristic spin-orbital excita- 
tions of finite width within the spin-orbital continuum. 
As shown in Fig.[TT]we can extract the dispersion of these 
characteristic excitations and find that the coupled spin- 
orbital excitations are gapless for the parameters con- 
sidered here. However, as can be clearly seen in Figs. 
E^b) andjSJd), the weights of the low-energy excitations 
are orders of magnitude smaller than the excitations lo- 
cated at higher energies. Hence the excitations at high 
energies give the most dominant contribution to the dy- 
namical spin structure factor. We therefore expect that 
these excitations will generate additional entropy in the 
corresponding temperature range which should show up, 




FIG. 12: (Color online) Dynamical spin structure factor 
S(q, ui) calculated perturbatively for the spin-orbital chain at 
temperature T/J = 0.1 with y = — 1. 



for example, in the specific heat which will be studied in 
the next section. 

Moreover, we find that these coupled spin-orbtial exci- 
tations are responsible for the Kohn anomalies mentioned 
above. We observe that the Kohn anomalies at interme- 
diate momenta occur when the energy of the spin wave 
coincides with that of a characteristic spin-orbital exci- 
tation. This is different from the Kohn anomaly in the 
spin-wave dispersion of itinerant ferromagncts. In this 
case the interaction between the localized spins given by 
the lattice ions is induced by scattering with conduction 
electrons and hence the Kohn anomaly is determined by 
the shape of the Fermi surface^ - — The Kohn anomaly 
we find within the present context is also due to interac- 
tion effects, where the nature of the interactions - coupled 
spin-orbital degrees of freedom - is distinct from the ones 
of the itinerant ferromagncts. For the case x = —y = 1 
the spin- wave dispersion has a discontinuity of the order 
Aw ~ 0.01J at the point q ~ 0.509tt (see inset of Fig.[SJ). 
However, for the crossing points of the magnon and the 
coupled spin-orbital excitation located at q ~ 0.377T and 
q ~ 0.76-7T (see Fig. [TTj) no Kohn-anomaly could be re- 
solved. We believe that this is a consequence of the 
weight of the coupled spin-orbital excitations: Whereas 
at q ~ 0.5097T the imaginary part of the self-energy dis- 
plays a steep increase of several magnitudes, at the other 
crossing points the slope towards the local maxima is far 
more moderate. 

At finite temperatures two effects contribute to the 
broadening of the central peaks of the dynamical spin 
structure factor. First, there is a broadening due to ther- 
mally excited magnons which is already present in the 
ID Heisenberg chain discussed in Sec. |TVl This is com- 
bined with the broadening due to the interaction with 
the orbital degrees of freedom. Here the BF continuum 
is smeared out by thermal fluctuations compared to the 
zero temperature case. Results for the structure factor 
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FIG. 13: (Color online) Magnon linewidth T q (FWHM of 
S(q,u))) as a function of q, as obtained at y = — 1 and tem- 
peratures T/J — 0, 0.1 and 0.2. The lines are guides to the 
eye. 

at finite temperatures are shown in Fig. 1121 

The broadening at small momenta is dominated by 
thermal fluctuations and the lineshape is very similar to 
that of the pure spin model discussed in Sec. HVl A 
further strong broadening in going from q = 0.37T to 
q = 0.57T signals the relevance of coupled spin-orbital 
degrees of freedom on the spin dynamics at intermedi- 
ate and high momenta. Again, additional structures in 
S(q, ui) are visible related to the spin-orbital excitations 
discussed above. 

Finally, we analyze the variation of the FWHM with 
increasing temperature for a representative value of y = 
— 1, sec Fig. [T3] At T > the thermal broadening at 
small momenta is clearly visible. As in the zero tempera- 
ture case another strong increase of T q between q = 2ir/5 
and q = ir/2 is observed due to coupled spin-orbital de- 
grees of freedom. For T/J~ 0.1 and T/J = 0.2 we 
observe that T q has a temperature dependent maximum 
from where T q decreases towards the boundary of the 
Brillouin zone. This is due to the fact that the ther- 
mal broadening of the central peaks of the dynamic spin 
structure factor decreases from intermediate to high mo- 
menta (see Fig.EJ). Actually without coupling to any or- 
bital degrees of freedom we expect r g to be very small at 
the boundary of the Brillouin zone. Thus a large band- 
width at q = 7r makes the spin-orbital model distinct 
from a pure ID Heisenberg ferromagnet. 



VI. THERMODYNAMICS 

Coupled spin-orbital degrees of freedom will not only 
influence the spin dynamics but also the thermodynam- 
ics of the system. We expect, in particular, that the 
spin-orbital excitations which were shown to affect the 




FIG. 14: Diagramatic representations of the second order con- 
tributions to the free energy as given by Eq. Q6.6JI . 

dynamical spin structure factor in the previous section 
will also become observable in thermodynamic quanti- 
ties when comparing the MF decoupling and the per- 
turbative solution. In order to investigate this issue we 
rewrite Eq. (|2.1j) as 

H St (1) = H mf + 5H, (6.1) 

with SH = H S t{1) - H M f- The MF part reads 

N 

H M F = Y] {Jt [Tj-T j + +J S Sj ■ Sj+1 

■j=i (6-2) 
- (Sj ■ S j+1 ) MF ([Tj ■ T j+1 ] 1 ) MF } . 

The exchange constants Js,t arc defined in Eqs. (|3.3|) . 67 
We use the Hamiltonian (|6.1jl to determine the free en- 
ergy per site perturbatively, following the expansion, 

/ = /mf + /mf + ( sh )'mf - ( Sh2 )mf + ■ ■ • ' 

(6.3) 

where /j^ F (/mf) ^ s the expression for the free energy per 
site stemming from the spin (pseudospin) sector within 
the MF decoupling solution. The subscript indicates that 
the respective correlation functions are calculated with 
Huf- Moreover the superscript c means that the above 
expansion of the free energy is restricted to connected 
diagrams. We note that Eq. (|6.3p is a high temperature 
expansion valid if \x\T/J > 1 and \y\T/J > 1. 

A straightforward calculation shows that the first or- 
der contribution only shifts the free energy, Eq. (|6 .3(1 . 
and will not show up in thermodynamic observables ob- 
tained by taking derivatives of the free energy. For the 
second order contribution we have to evaluate two- and 
four-point correlation functions both for the spin and the 
pseudospin part. We use the abbreviations 

((S r S j+1 )(S r S l+1 )) = a + b(j,l) (6.4) 

for the spins and 

([Tj ■ T j+1 ] x [ti -ti+xIj) =c+d(j,l) (6.5) 

for the pseudospins. Here the site-independent quanti- 
ties a and c stand for the disconnected parts of the four- 
point correlation functions whereas b(j, I) and d(j, I) fol- 
low from the connected ones. One finds after a straight- 
forward calculation that only the product of the con- 
nected parts contributes to the second order correction, 
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leading to 



N 



(5H>) 



MF 



J 2 l)d(j, I) . (6.6) 

3,1=1 



To proceed further, we again apply the MSWT to the 
spin and a Jordan- Wigner transformation to the pseu- 
dospin part. The evaluation of d(J, I) is again straight- 
forward and yields 



d(j,l) 



X - 2 n F [ W r(fci,0)]{l-n F [ W r(fe,0)]} 



2N 2 



ki,k 2 



xe i(ki-k 2 )(j-l) {i + co ^ fci + j^y 



(6.7) 

with ujp F (q, 5) as given in Eq. (|3T6j) . 

The evaluation of Eq. f)6.4|) is more involved. We first 
apply a Dyson-Maleev transformation and treat the ob- 
tained expressions using Wick's theorem. In addition, we 
also have to account for the constraint of nonzero mag- 
netization at T > imposed by the MSWT. The cor- 
responding diagrams are shown in Fig. 1141 Within this 
approximation the specific heat per site reads 



with 



c = c M f + c 2 + 



N 



C2 



r, d 2 j,l=l 



Q T 2 



2NT 



(6.9) 



We calculate the first term in Eq. (|6.8[) within the MF 
decoupling, cmf = c^p + c mfj from the internal energy 
which is determined by the respective nearest-neighbor 
correlation functions allowing us to keep terms up to 
quartic order in the bosonic operators^ This strategy 
makes it possible to obtain reliable results for cfj F up to 
T/(| Js\S 2 ) < 1. The second order correction c 2 given in 
Eq. (|6.9p is obtained using the Dyson-Maleev transfor- 
mation so that quartic terms are also included and the 
order of approximation is the same. Since we are using 
a high temperature expansion, Eq. (|6.3p , in combination 
with the MSWT to evaluate the diagrams, our results 
are only valid in an intermediate temperature regime. If 
we restrict ourselves to parameters x — — y > then this 
temperature range is given by 1/x <C T/J <C xS 2 . In 
the following we therefore only consider the case x 3> 1 
and compare the results from perturbation theory with 
numerical data obtained by TMRG. 

As shown in Fig. [121 the specific heat c/(Jx) 2 exhibits 
a broad maximum which corresponds to the character- 
istic energies of spin and fermionic particle-hole excita- 
tions. In the temperature range where the perturbativc 
approach is valid we find excellent agreement with the 
numerical solution. In particular, the perturbative cor- 
rection C2 (see Fig. [IB")) correctly captures the weight shift 
from low to intermediate temperatures visible when com- 
paring the numerical and the MF decoupling solution. In 




0.4 0.6 
T/(Jx) 



FIG. 15: Specific heat per site c/(Jx) 2 as a function of T/Jx 
for (a) x = 10 and (b) x — 2. The circles denote the nu- 
merical data from TMRG with the solid line obtained by 
a low-temperature fit of the TMRG data for the inner en- 
ergy. The dashed lines correspond to MF decoupling while 
the dashed-dotted lines are the results obtained by perturba- 
tion theory. The perturbative results are expected to be valid 
for 1/x 2 <C T/(Jx) <C 1. Insets: Specific heat cmf within the 
MF solution (dashed-dotted line) with the contributions from 
the spin (cmf solid line) and the orbital (cjjp dashed line) 
sector shown separately. 



spite of this weight shift, the MF decoupling yields overall 
a very reasonable description of the specific heat for both 
cases shown in Fig. [TSJ For x = 10, (Fig. [T5ja)) the spe- 
cific heat c has a broad maximum at T / (Jx) ~ 0.4. This 
maximum results from a distinct maximum in the orbital 
contribution cJj F , see insets in Fig. 1151 In contrast the 
spin contribution cj^ F increases steadily with increasing 
temperature, in agreement with the higher energy scale 
for spin excitations. As a result, the total specific heat 
has only a weaker and broader maximum than suggested 
by the orbital part. 

The MF decoupling does seem to fail, however, at very 
low temperatures. Here the MF solution predicts that 
spin excitations give the dominant contribution leading 



15 



X 
— i 

-0. 



1 1 1 1 1 1 1 1 1 

' ^ 

/ \ 
// \ 




/ 1 




/ / 




/ ; 




/ / 




/ / 


— 


/ ' 




/ / 




/ / 




/ 1 




I i 1 1 1 





0.2 0.4 0.6 0.8 



T/(Jx) 

FIG. 16: (Color online) Perturbative contribution to the spe- 
cific heat per site C2/(Jx) 2 , Eq. (|6.9p . as a function of T/Jx 
for x — 10 (solid line) and x = 2 (dashed line). 



to a c(T) ~ yf behavior. This is in contrast to an ex- 
trapolation of the numerical data shown as dot-dashed 
lines in Fig. fTS] which suggests an approximately linear 
dependence on temperature. This discrepancy comes as a 
surprise because our perturbative calculations of the dy- 
namical spin-structure factor in Sec. IV Bl lead us to the 
conclusion that the magnons survive as sharp quasipar- 
ticles at low energies. More generally, one might argue 
that the ground state does not show any entanglement 
between the two sectors because of the classical nature 
of the FM state thus allowing for spin-wave excitations 
at low energies. From the point of view of a low-energy 
effective field theory, however, the situation is much less 
clear. While a FM chain is described by a low-energy 
effective theory with dynamical critical exponent 2 = 2, 
the fermionic orbital chain has z = 1. A coupling of 
spatial spin deviations to time-dependent orbital fluc- 
tuations then seems to require that the low-energy ef- 
fective theory for the coupled system has 2 = 1. As a 
consequence the temperature dependence of c(T) would 
indeed be linear. However, such an approach leaves open 
the role and treatment of the Berry phase terms. 

Within the perturbative approach the open question 
is whether or not higher order contributions to the self- 
energy might induce a significant broadening of the dy- 
namical spin-structure factor also at low energies. In this 
regard we note that the vertex responsible for the broad- 
ening of S(q, u) studied in Sec. IV Bl docs not play any 
role for the thermodynamics of the system. Here the con- 
straint of vanishing magnetization means that such dia- 
grams do not contribute to static correlation functions so 
that the lowest order corrections are caused by the vertex 
shown in Fig. [H] 



VII. CONCLUSIONS 

In summary, we have investigated coupled spin-orbital 
degrees of freedom in a one-dimensional model. For fer- 
romagnetic exchange we have shown that the considered 
model at special points in parameter space can be written 
in terms of Dirac exchange operators for spin S and pseu- 
dospin t. As a consequence, three collective excitations 
of spin, orbital and coupled spin-orbital type do exist. In 
particular, we discussed the case of Dirac exchange oper- 
ators for S = t = 1/2 where the dispersions of all three 
elementary excitations are degenerate and lie within the 
spin-orbital continuum. While the spin-orbital excita- 
tions stay confined in this decay becomes possible 
once we move away from this special point. 

For antifcrromagnetic exchange the one-dimensional 
spin-orbital model captures fundamental aspects of 
physics relevant for transition metal oxides and, as we 
have shown, sharp excitations do not exist. To address 
the question how the spin dynamics is influenced by fluc- 
tuating orbitals we considered the extreme quantum limit 
of orbitals interacting via an XY-type coupling. This al- 
lowed us to map the orbital sector onto free fermions 
using the Jordan- Wigncr transformation. In spin-wave 
theory the spin sector is described by bosons so that our 
model corresponds to an effective boson-fermion model 
which applies for low temperatures. An analytic calcu- 
lation of the properties within a mean-field decoupling 
approach is then straightforward. Compared to a numer- 
ical phase diagram based on the density-matrix renor- 
malization group we find that the regime with ferromag- 
netically polarized spins is stabilized by the decoupling 
procedure. Furthermore, the mean-field decoupling gives 
rise to a finite temperature dimerized phase for certain 
parameters when starting from the ferromagnetic ground 
state. While a phase with long-range dimer order at finite 
temperatures is not possible in a purely one-dimensional 
model, the mean-field approach also completely ignores 
any kind of coupled spin-orbital excitations. 

Thus we developed a self-consistent perturbative 
scheme to explore the role played by spin-orbital cou- 
pling. In perturbation theory the boson-fermion inter- 
action does not produce any infrared divergencies in the 
one-dimensional model due to the lack of nesting. This 
makes a perturbative calculation of the spin structure 
factor S(q, tu) possible. At large momenta q, we find that 
S(q, w) shows a significant broadening due to scattering 
of magnons by orbital excitations. For small momenta, 
on the other hand, no broadening in this lowest order 
perturbative approach is observed because the magnon 
cannot scatter on these excitations. The onset of the 
broadening occurs at momenta where the magnon enters 
the boson-fermion spectrum. This point, as well as de- 
tails of the full width at half maximum is determined by 
the strength of interaction. Most interestingly, S(q, u>) 
does show additional peaks and shoulders corresponding 
to characteristic spin-orbital excitations. At points where 
the renormalized spin- wave dispersion and the dispersion 
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of these excitations cross, Kohn anomalies do occur. 

Furthermore, we compared numerical data for the spe- 
cific heat of the spin-orbital model with the mean-field 
decoupling solution and an approach where we also took 
the second order correction to the mean-field result into 
account. Overall, we found that the mean- field decou- 
pling does describe the specific heat reasonably well. A 
redistribution of entropic weight from low to intermedi- 
ate temperatures observed when comparing the numeri- 
cal data and the mean-field solution is very well captured 
by the second order perturbative correction. An interest- 
ing open point is the behavior of the specific heat c(T) at 
low temperatures. While the mean-field solution predicts 
c(T) ~ VT due to spin-wave excitations, the numerical 
data suggest instead that c(T) ~ T. We have speculated 
that a coupling of the two sectors might indeed lead to a 
low-energy effective theory with dynamical critical expo- 
nent z = 1 but details of such a theory need to be worked 
out in the future. 

In conclusion, we have shown that while collective 
spin-orbital excitations with infinite lifetime do not ex- 
ist for antiferromagnetic superexchange the coupled spin- 
orbital degrees of freedom have a strong influence on the 
spin excitation spectrum as well as on the thermody- 
namic properties of the system. However, the treatment 
of spin-orbital systems beyond the range of validity of 
mean-field decoupling and perturbative schemes remains 
an open problem in theory. 
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Appendix: Details of the perturbative approach for 
the boson-fermion model 



and 



X +Jj - cos fe)^b,fc j cos q. (A. 2) 

k / 



w f(<7) = yS 2 + x ■ 

For T = only the magnon dispersion is renormalized to 

C(9) = 2JS(|y| + l/7r)(l-co S(? ). 

By virtue of Dyson's equation we may calculate the 
bosonic Green's function at T = Since the MF de- 
coupling already takes the first order contributions into 
account, the lowest order diagrams we obtain are of sec- 
ond order. The self energy is thus approximated by 
the proper self-energy obtained by summing up those 
diagrams which may be composed by the second order 
diagrams. From this we have £(q,w) w £2(9,0;) = 
Sg F (g,a;) + Eg F (<7,aj) + £2,1(9) where the diagrams are 
given by £+ p (g,w) (Fig.Efa)), E^ F (q,cj) (Fig.^b)), and 
£2,1(9) (Fig.[6fc)). A straigthforward calculation reveals 
that at zero temperature £2,1(9) and £g P (<jr, w) vanish. 
For £g F (q,oj) we find the expression given in Eq. (|5.5|) . 

For finite temperatures we calculate the imaginary 
time Green's function. We find 

r 



a . b 



^G F \k2,UJa,F)G F 0) (q — ki + fc 2 ,W&, F ) 

(A.3) 



X G^ (h,±(Uv,B ~^b,F +W ttjF )) 



for the diagrams in Fig. [5Ia,b) , and 



T 

£2,1(9) = — ttj 2^ u BV (q 1 q,k 1 )u BV (k 2l k2ik 1 ) 



N 2 



[G^ih,^) 



2 ( A - 4 ) 



We start by a MF decoupling, rewriting the interaction 



H 1 = H 



1,MF 



(#1 - Hi. 



MF 



SH 



with 



#1,MF ^jp^2^BF(k,k,q) 2 n b , k f%f q + nf,qb],b k 



We treat SH as perturbation. The averages n aiP = 
(ata p ) are determined self-consistently within this MF 
scheme. This leads to a renormalization of the magnon 
and fcrmion dispersions. We find 

C(<?) =2JS ^\y\ - j^Yl cosk C^-^l) ~ A*) 

(A.l) 



for the diagram given in Fig. [Sic). Here we have used 
5 F 0) (9,w MiF ) = (iu^ F - LUFiq))- 1 and G^\q,^ v .B) = 
{iuj u ,B — C(9)) _1 as the fermionic and bosonic Matsub- 
ara Green's function for the nonintcracting Hamiltonian, 
respectively. w„. F are the fermionic Matsubara frequen- 
cies. After performing the frequency sums we end up 
with 



s 2,i(g) 



1 



TN 2 ^ 

k\ ,k 2 



wbf(<7, q, h)uiBF(k2,k2, fci) 



x n B [C(fc2)](l - n F [w F (fci)])n F [w F (£:i)], 



(A.5) 

and the self-energies Ti BF (q,u} Vt B) as given in Eqs. (|5.7p 
and (JlU). 
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